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1.  Introduction 


Three-dimensional  (3-D)  textile  composites  offer  key  benefits  to  structural  robustness  for 
applications  involving  impact,  multidimensional  loading  or  thick  sections  with  relatively  large 
through-thickness  or  delamination  stresses.  Such  is  also  the  case  for  any  two-dimensional  (2-D) 
woven  or  braided  composite  structures,  which  inherently  provide  reinforcement  in  multiple 
directions  and  may  include  some  out-of-plane  reinforcement  by  nature  of  the  undulation  of 
interwoven  fiber  tows.  However,  the  3-D  orthogonal  weave  incorporates  fiber  tows  directly  in 
the  through-thickness  direction.  This  is  critical  to  delamination  resistance,  an  otherwise  common 
weak  point  of  thick  composite  structures.  Consequently,  this  also  leads  to  loss  of  in-plane 
properties  compared  to  a  2-D  or  laminated  structure  of  the  same  fiber  volume  fraction,  as  the 
stitches  lead  to  interruption  of  in-plane  fiber  tows  and  loss  of  volume  fraction  of  fibers  aligned  to 
resist  in-plane  loading. 

The  increased  microstructural  complexity  of  textile  composites  also  leads  to  increased 
complexity  of  characterization  and  analysis  as  well  as  a  need  for  nontraditional  analysis  methods. 
The  benefits  in  performance  and  manufacturability  can  outweigh  these  drawbacks  with  effective 
design  and  knowledge  of  the  micromechanical  response,  as  appropriate  to  the  particular 
application  being  considered. 

An  accurate  model  must  accommodate  the  fact  that  3-D  woven  composites  exhibit  multiple 
potential  failure  mechanisms  (7),  which  depend  upon  the  loading  conditions  and  particulars  of 
the  layup  and  materials.  Three-dimensional  weaves  consistently  show  improved  damage 
resistance  over  their  2-D  counterparts  due  to  the  energy  absorbing  capacity  enhanced  by  the 
z-direction  fiber  tows  (2,  2).  Three-dimensional  weaves  also  show  more  capacity  to  absorb 
multiple  strikes  before  perforation  and  show  less  damage  localization  ( 4 ).  This  can  be  of  particular 
interest  in  armor  applications,  and  related  modeling  has  addressed  issues  surrounding  the 
computational  needs  required  to  reflect  the  energy  absorption  modes  and  damage  progression  (5). 

In  addition  to  the  improvement  of  impact  performance,  the  effects  of  z-stitching  upon  shear 
properties  have  been  shown  experimentally,  with  increased  delamination  resistance  and 
significantly  increased  compression-after-impact  (CAI)  capacity  ( 6 ,  7).  Failure  modes  are  shown 
to  vary,  dependent  on  the  particulars  of  the  textile  architecture,  especially  with  regards  to  CAI. 

Designs  of  current  textile  composite  structures  are  often  based  upon  well-known 
phenomenological  failure  criterion,  predominantly  the  maximum  stress  criterion,  maximum 
strain  criterion,  and  quadratic  interaction  criterion,  such  as  the  Tsai-Hill  and  Tsai-Wu  failure 
theories  ( 8 ,  9).  Dedicated  analysis  of  2-D  and  3-D  textile  composites  has  been  approached 
before  (10-16);  however,  effective  strength  prediction  is  still  a  topic  of  much  current  study  and 
development. 
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Finite-element  model  (FEM)  micromechanical  methods  for  strength  modeling  of  textile 
composites  have  been  explored  in  previous  works  (1 7-19).  This  employed  a  method,  known  as 
the  Direct  Micromechanics  Method  (20,  21),  to  develop  failure  envelopes  and  a  quadratic  stress 
gradient  failure  theory  for  a  plain-weave  textile  composite.  Therein,  it  was  shown  that  the 
microlevel  scale  of  a  textile  composite  is  physically  larger  and  geometrically  more  complex  than 
a  unidirectional  composite.  Thus,  many  common  assumptions  and  traditional  micromechanical 
analysis  techniques  break  down.  Improvements  to  these  techniques  have  been  presented  towards 
effective  failure  prediction  under  complex  loading  conditions. 

Three-dimensional  woven  composites  show  solid  potential  for  use  in  improving  delamination 
resistance  and  impact  resistance  in  a  lightweight  material,  which  can  be  important  in  armor  or 
structural  applications.  In  particular,  with  the  electromagnetic  railgun  launcher,  significant 
electromagnetic  forces  generated  in  the  aluminum  rail  sections  can  lead  to  delamination  of  the 
insulator  materials,  currently  a  filament-wound  S2  glass  -  BMI  composite.  Therefore, 
investigation  is  made  into  the  delamination  resistance  potential  of  composites  with  3-D 
reinforcement.  The  current  work  involves  applying  FEM  micromechanics  to  the  3-D  orthogonal 
weave  to  provide  detailed  stiffness  and  strength  modeling  for  failure  prediction.  Experimental 
analysis  of  multiple  specimen  types  has  been  performed  to  provide  model  verification.  The 
established  micromechanical  model  is  then  applied  towards  a  parametric  study  of  the  effects  of 
variation  in  stiffness  and  strength  for  various  stitch  densities,  which  is  used  to  illustrate  the 
material’s  potential  and  provide  guidelines  for  microstructural  optimization. 


2.  Methods 


Finite-element  (FE)  based  micromechanical  analysis  has  been  performed  to  determine  the 
stiffness  and  strength  of  a  3-D  orthogonal  weave  textile  composite.  Once  characterization 
simulations  have  been  performed  within  the  material  model,  FE  stress-field  outputs  become  the 
basis  for  an  analytical  model  which  can  predict  failure  under  arbitrary  loading  conditions.  This 
then  allows  for  a  parametric  study  in  which  the  properties  of  various  3-D  orthogonal  fiber 
microarchitectures  can  be  investigated  towards  determination  of  an  optimal  geometry. 

2.1  Representative  Volume  Element  (RVE)  Development 

In  terms  of  structural  layout  and  dimensional  accuracy,  an  extensive  series  of  photomicrographs 
was  taken  to  ensure  the  precise  representation  of  the  microstructure  in  the  FEM.  Specimens 
were  cut  in  several  planes  of  the  3-D  orthogonal  source  plaque  at  various  magnifications,  with  a 
length  scale  attached  by  the  photomicrograph  software.  Figure  1  shows  a  picture  taken  in  the  x-y 
plane  (a  top  view)  from  which  multiple  dimensions  can  be  ascertained,  such  as  stitch  cross-section, 
stitch  spacing,  in-plane  tow  width,  etc.  The  photo  of  figure  2  is  taken  in  th ey-z  plane  (a  side  view), 
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Figure  1 .  Photomicrograph  in  the  x-y  plane. 


Figure  2.  Photomicrograph  in  the  y-z  plane. 
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which  illustrates  the  warp  and  weft  tow  heights,  as  well  as  multiple  matrix  pockets  in  areas 
where  tows  cannot  be  strung  or  woven  (in  this  case,  due  to  the  interference  of  stitching).  Further 
pictures  not  shown  here  were  also  taken  to  complete  the  microarchitectural  and  dimensional 
determination  as  well  as  to  provide  assurance  of  consistency  of  dimensions  at  multiple  locations 
along  the  material  plaque.  Figure  3  provides  a  postbumout  photomicrograph  for  which  all  resin 
material  has  been  removed  by  residence  in  a  high-temperature  furnace.  This  is  done  to  ensure 
proper  distinction  of  the  tow  and  matrix  phases  as  well  as  to  provide  an  accurate  measure  of  the 
total  fiber  volume  fraction.  These  pictures  are  not  used  to  determine  dimensions,  as  some  drift  of 
the  fiber  tows  will  occur  while  not  contained  in  the  binding  resin  matrix  (some  bowing  of  the 
tows  can  be  seen  in  figure  3,  which  does  not  occur  in  the  cured,  resin-impregnated  structure). 


Figure  3.  Photomicrograph  in  the  x-y  plane  after  resin  matrix  burnout. 


Figure  4  shows  one  of  several  photomicrographs  taken  at  higher  magnification  to  investigate 
some  of  the  subtleties  of  the  microarchitecture.  Fiber  kinking,  which  is  sometimes  a  concern  for 
textile  composites,  is  seen  to  be  minimal  or  negligible.  Further,  even  at  this  high  magnification, 
the  orthogonality  of  tows  is  quite  evident.  In  fact,  there  is  surprisingly  little  “blunting”  or 
rounding  at  the  comer  points  at  which  tows  press  against  neighboring  tows.  Whereas,  for  textile 
composites,  in  general,  loss  of  ideal  geometry  can  sometimes  be  seen  in  the  postprocessing  (or 
postmanufacture)  fiber  structure,  especially  when  the  fiber  preform  is  draped  over  a  part  or  mold 
(though,  in  this  case,  the  source  material  is  a  straight  plaque). 
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Figure  4.  High  magnification  photomicrograph  illustrating  orthogonality  of  tows. 


The  smallest  repeatable  RVE  of  the  textile  architecture  under  investigation  is  shown  in  figure  5, 
with  dimensions  listed  in  table  1 .  Constituent  properties  are  shown  in  tables  2  and  3  (note  that 
fiber  tows  have  the  properties  of  an  orthogonal  composite  at  the  micro  or  element  level).  In 
figure  5,  the  “warp”  tows  represent  fiber  tows  for  which  the  fibers  run  in  the  direction  of  the 
x-axis.  “Weft”  tows  represent  those  for  which  fibers  run  in  the  y-direction.  Unlike 
microgeometry  such  as  the  2-D  plain  weave,  there  is  no  undulation  or  interweaving  of  the  fiber 
tows  (though  such  3-D  architectures  do  exist).  Thus,  each  contributes  its  strength  fully  to  the 
in-plane  properties  of  the  composite,  while  z-direction  tows  or  stitches  are  placed  to  provide 
through-thickness  properties.  During  fabrication  of  the  3-D  orthogonal  preform,  small  gaps 
subject  to  the  limitations  of  the  inter-strung  fiber  tow  exist.  Thus,  some  areas  of  interstitial 
matrix  will  be  present  (also  indicated  in  figure  5).  Given  these  interstitial  gaps  as  well  as  the 
inherent  volume  fraction  of  resin-impregnated  fiber  tows,  the  overall  volume  fraction  of  the 
current  RVE  is  54%. 

2.2  Modeling  for  Stiffness  Prediction 

The  RVE  is  subjected  to  individual  macroscopic  (or  RVE-level)  displacements  utilizing  periodic 
boundary  conditions  which  correspond  to  unit-strain  along  each  axis.  Each  independently 
applied  unit-strain  results  in  a  complete  microlevel  (or  element-level)  stress  field  comprised  of 


Figure  5.  RYE  geometry. 


Tablet.  RYE  dimensions. 


Dimension 

(mm) 

Description 

Value 

(mm) 

a 

RVE  width 

7.25 

b 

RVE  depth 

6.0 

h 

RVE  height 

1.68 

h 

Warp  (half)  thickness 

0.49 

h 

Weft  thickness 

0.70 

Wj 

Warp  width 

6.62 

w2 

Weft  width 

4.75 

Table  2.  Fiber  tow  and  matrix  material  properties. 


Material 

Ex 

(GPa) 

E2 

(GPa) 

Gn 

(GPa) 

Vl2 

S2  glass  composite 
(65%  fiber  volume) 

71.1 

16.2 

6.46 

0.25 

Neat  resin 

4.60 

4.60 

1.76 

0.31 

Table  3.  Fiber  tow  and  matrix  strength  properties. 


Material 

sL<+) 

(MPa) 

sL(-> 

(MPa) 

c  (+) 

(MPa) 

Sth 

(MPa) 

Slt 

(MPa) 

S2  glass  composite 

1705 

550 

45 

150 

50 

Neat  resin 

103 

103 

103 

103 

— 
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the  3-D  stresses  in  each  element.  The  RVE  is  then  characterized  by  analysis  of  the  macrolevel 
stresses  across  the  entire  RVE,  which  are  calculated  by  weighted  averaging  of  the  microstresses 
across  the  RVE  volume.  By  applying  known  independent  unit  strain  cases  and  evaluating  the 
resulting  macrostresses,  one  can  completely  determine  the  constitutive  matrix.  In  this  approach, 
structural  stiffness  coefficients  are  computed  directly  from  the  micromechanical  models,  and  no 
underlying  assumptions  or  approximations  are  needed. 

The  RVE  is  subjected  to  macroscopic  stresses  (a)  which  are  related  to  macroscopic  strains  (e) 
according  to  the  following  equation: 

[<t]  =  [C]{£}.  (1) 

The  constitutive  matrix  [C]  must  be  evaluated  to  define  this  correlation.  To  completely 
determine  [C],  a  single  unit  strain  is  applied  and  all  other  deformation  terms  are  set  to  zero  in 
each  of  the  cases  shown  in  table  4.  Displacements  are  applied  using  periodic  boundary 
conditions.  The  periodic  displacement  boundary  conditions  isolate  the  mechanical  effects  of 
application  of  unit  strains  or  curvatures  and  ensure  the  repeatability  of  deformations.  Thus,  the 
RVE  is  not  only  repeatable  as  a  representative  geometry,  but  is  also  mechanically  repeatable  in 
that  each  RVE  has  an  identical  response  to  strains  and  curvatures  regardless  of  the  location  of 
that  RVE  in  a  textile  plate  or  component.  Note  that  displacements  u,  v,  or  w  correspond  to  axial 
directions  x,  y,  and  z,  respectively. 


Table  4.  Periodic  displacement  boundary  conditions. 


u(b,y,z)- 

u{0,y,z) 

v(b,y,z)- 

v(0,y,z) 

w(x,y,h)- 

w(x,y,0) 

u(x,a,z)- 

u(x,0,z) 

v(x,a,z)- 

v(x,0,z) 

w(x,a,z)- 

w(x,0,z) 

w(x,b,z)- 

w(x,0,z) 

£*=  1 

b 

0 

0 

0 

0 

0 

0 

£y=  1 

0 

0 

0 

0 

a 

0 

0 

£z=  1 

0 

0 

h 

0 

0 

0 

0 

Jxy=  1 

0 

a/2 

0 

b/2 

0 

0 

0 

Y*z=  1 

0 

0 

0 

0 

0 

b 

0 

V=1 

0 

0 

0 

0 

0 

0 

a 

The  FEM  results  for  each  element  yield  the  microstresses  resulting  from  an  applied  macrolevel 
strain.  The  corresponding  macrolevel  stress  in  each  case  can  be  computed  by  averaging  the 
microstresses  over  the  entire  volume  of  the  RVE  as  follows: 

(2) 

Here,  e  denotes  summation  over  all  elements  in  the  FE  model  of  the  RVE,  Vs  is  the  volume  of 
each  element  e,  and  a,  b,  and  h  are  the  dimensions  of  the  RVE. 
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Thus,  the  constitutive  relation  of  equation  1  can  be  found  by  independently  evaluating  the  strain 
cases  shown  in  table  4  and  incorporating  FEM  results  per  equation  2.  Each  case  will  determine 
one  column  of  the  stiffness  matrix  [C],  which  is  then  completely  defined  once  the  results  of  all 
cases  have  been  combined. 

In  order  to  then  evaluate  the  stress  state  in  the  RVE  resulting  from  any  general  applied  load  case, 
no  further  FEM  analysis  is  necessary.  The  microstresses  in  each  element  can  be  extrapolated 
from  the  preliminary  RVE  analysis  (just  described)  of  each  of  the  linearly  independent 
macrostrain  components.  The  microstress  state  for  a  general  applied  stress  is  obtained  by 
superposing  multiples  of  the  results  from  the  unit  macrostrain  analysis  as  follows: 

{<r*}  =  [F*]{c}  .  (3) 

Here,  the  6x6  matrix  [E6]  contains  the  microstress  in  each  element  resulting  from  the  unit  strain 
analysis  cases  of  table  4.  For  example,  the  microstress  <jy  in  the  RVE  for  £x  =  0.05  and  £y  =  0.003 
is  calculated  as  <7y  =  O.O5F21  +  O.OO3E25.  Thus,  further  stress  analyses  are  performed  in  an 
analytical  model  based  upon  superposition  of  the  results  of  unit  strain  FEM  analyses. 

2.3  Strength  Prediction 

Strength  can  then  be  determined  by  evaluating  the  microstress  field  in  the  RVE  resulting  from 
any  prescribed  macrolevel  force  or  stress.  Strength  is  predicted  by  comparing  the  computed 
microstresses  in  each  element  against  failure  criteria  for  the  constituent  materials.  Failure  is 
checked  on  an  element-by-element  basis,  and  the  failure  criterion  of  each  element  can  be 
selected  appropriately  based  upon  whether  it  is  a  tow  or  matrix  element.  For  the  isotropic  matrix 
elements,  the  maximum  principal  stress  criterion  has  been  chosen  to  evaluate  element  failure. 

For  fiber  tow  elements,  the  Tsai-Wu  failure  criterion  is  used,  as  this  criterion  is  more  suited  to 
the  orthotropic  nature  of  the  tows  including  stress  interaction  effects. 

Failure  envelopes  (or  simpler  case-by-case  failure  analyses  appropriate  to  an  expected  load  set) 
are  generated  by  first  selecting  a  macrostress  state  to  investigate.  Then,  the  macrostrains  and 
curvatures  resulting  from  this  applied  loading  are  calculated  from  equation  1 .  Based  on  the 
scaled  superposition  of  the  results  from  FEM  analysis  of  the  unit  load  cases  shown  in  table  4,  the 
resulting  stress  field  for  the  entire  RVE  is  then  calculated  by  equation  3.  Failure  is  then  checked 
in  each  element  against  a  given  failure  criterion.  This  cycle  is  then  repeated  while  progressively 
increasing  a  selected  force  or  moment  resultant  and  holding  all  others  constant  until  an  element 
level  failure  criterion  is  exceeded.  If  a  particular  failure  criterion  is  exceeded,  the  element  and 
the  RVE  are  considered  failures,  which  then  defines  the  threshold  of  the  failure  envelope  at  a 
given  point.  Thus,  failure  envelopes  for  the  textile  composite  can  be  generated  in  various  force 
and  moment  resultant  spaces. 
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2.4  Experimental  Methods 

Several  experimental  tests  have  been  performed  to  verify  the  accuracy  of  stiffness  and  strength 
predictions  via  the  micromechanics  modeling.  This  establishes  a  baseline  of  veracity  for  the 
modeling  and  also  establishes  the  material  behavior  in  terms  of  linearity  vs.  nonlinearity, 
catastrophic  vs.  progressive  failure,  failure  modes,  etc.  Two  test  sets  have  been  performed — 
standard  tensile  tests  to  determine  stiffness  and  strength  in  the  axial  directions  and  Iosipescu 
shear  tests  to  determine  in-plane  and  transverse  shear  stiffness  and  strength.  To  ensure 
repeatability  and  statistical  certainty,  each  test  is  run  for  eight  specimens.  Tests  are  performed 
on  a  hydraulic  mechanical  testing  system  machine  with  a  crosshead  mounted  load  cell.  Full-field 
strain  measurements  are  obtained  via  the  digital  image  correlation  optical  measurement 
technique.  In  this  technique,  a  random  speckle  pattern  is  applied  to  the  gauge  section  of  the 
specimen.  A  pair  of  digital  cameras  then  records  a  series  of  stereo  images  to  track  changes  in  the 
speckle  pattern  during  testing.  The  images  are  postprocessed  with  image  correlation  software 
(22),  which  tracks  the  relative  displacement  of  all  speckles  within  the  pattern  and  computes  3-D 
surface  strains  from  these  displacements. 

In  addition  to  the  characterization  tests,  30°  and  45°  off-axis  specimens  have  been  tested.  These 
provide  more  challenging  test  cases  for  model  validation,  with  a  more  complicated  stress  state. 

As  shown  in  figure  6,  standard  uniaxial  specimens  are  machined  at  30°  or  45°  from  the  material 
axes  of  figure  5.  When  the  specimen  is  loaded  in  tension,  the  stress  state  seen  by  the  material 
(subscript  “ mcs ”  in  equations  4  and  5)  according  to  simple  tensor  change  of  axis  transformation 
(9)  indicates  that  axial,  transverse,  and  shear  stress  will  be  present. 


r 

0.5' 
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II 

0 
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^  mcs 

0.5  > 
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2.5  Parametric  Modeling  of  Tow  Microgeometry 

Once  the  FEM  micromechanical  tools  have  been  established  and  experimental  results  have 
validated  the  performance,  the  methods  can  then  be  extended  to  perform  a  parametric  study  of 
the  effects  of  changes  to  microgeometry  on  mechanical  performance.  This  becomes  a  powerful 
tool  by  which  the  effects  of  changes  to  microgeometry  on  the  mechanical  performance  of  the 
material  can  be  quickly  explored. 

The  chief  driving  force  behind  selection  of  a  3-D  weave  is  to  improve  through-thickness 
properties,  which  are  predominantly  affected  by  the  presence  of  the  tows  or  stitches  in  the 
thickness  direction.  Thus,  in  the  current  study,  the  effects  of  variation  in  stitch  density  (the 
center-to-center  distance  between  stitches,  as  indicated  by  the  number  of  stitches  per  centimeter) 
is  investigated,  in  terms  of  its  consequence  on  improvement  in  transverse  shear  stiffness  and 
strength,  as  well  as  through-thickness  properties.  Because  in-plane  properties  can  be  adversely 
affected  by  the  presence  of  stitching,  this  consequence  is  also  explored. 

Referring  again  to  figure  5,  the  parametric  study  is  performed  by  modifying  the  stitch  spacing. 

In  each  case,  a  new  RVE  is  created  to  reflect  the  new  stitch  density,  and  FE  micromechanics  are 
employed  to  determine  the  structural  response.  Increasing  the  RVE  width  and  depth  (dimensions 
a  and  b,  respectively)  corresponds  to  a  lower  stitch  density  and  vice  versa.  Tow  widths  w  are 
consequently  affected  as  well,  but  all  other  dimensions  remain  the  same.  Four  total  stitch 
densities  are  investigated,  as  shown  in  table  5.  Note  that  the  material  with  zero-stitch  density 
essentially  amounts  to  a  [0/90/0]s  cross-ply  laminate  since  only  in-plane  fibers  in  the  x  and  y 
axial  directions  are  present.  Because  the  properties  of  the  3-D  orthogonal  structure  can  be 
compared  to  those  of  a  cross-ply  laminate,  this  provides  a  further  basis  upon  which  results  may 
be  gauged. 
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Table  5.  Variation  in  RYE  geometry  for  parametric  study. 


Name 

Stitch  Density 
(stitches/cm) 

RVE  Width 
(mm) 

RVE  Depth 
(mm) 

Unstitched 

0 

9.0 

9.0 

Sparse 

1.1 

9.0 

9.0 

Baseline 

1.5 

7.25 

6.0 

Dense 

2.0 

5.0 

5.0 

Once  the  RVE  geometry  and  mesh  have  been  constructed,  the  methods  of  section  2. 1  can  then  be 
reapplied  for  each  new  RVE  in  table  5.  Displacements  corresponding  to  in-plane  and  or 
transverse  shear  can  be  applied  using  periodic  boundary  conditions.  The  resulting  stress  field  in 
each  case,  which  now  varies  as  the  microgeometrical  dimensions  change,  then  relates  the  known 
displacement  to  the  corresponding  altered  elastic  constant.  Also,  per  section  2.1,  strength  is 
again  determined  by  evaluating  the  stress  field  on  an  element-by-element  basis  to  determine  the 
maximum  applicable  macrolevel  stress. 


3.  Results 


3.1  Stiffness  and  Strength  Properties 

Stiffness  values  have  been  determined  from  the  modeling  and  experimental  methods  just 
described.  Although  modeling  can  quickly  predict  strength  under  any  loading  condition, 
logistical  limitations  dictate  that  experimental  tests  have  been  performed  to  determine  strength 
under  uniaxial  loading  conditions.  Though  not  shown  in  test-by-test  details,  each  of  the 
experimental  tests  from  table  6  has  been  repeated  for  eight  specimens.  An  average  value  is 
presented. 

Table  6.  Experimental  and  predicted  stiffness  and  strength. 


Description 

Experimental 

(Pa) 

Predicted 

(Pa) 

Difference3 

(%) 

Ex 

x-axis  stiffness 

2.25  elO 

2.41  elO 

7.3 

XT 

x-axis  strength 

5.48  e8 

5.65  e8 

3.1 

Ey 

y-axis  stiffness 

2.26  elO 

2.45  elO 

8.8 

Yt 

y-axis  strength 

4.31  e8 

4.96  e8 

15.0 

Gxy 

Shear  stiffness  (in-plane) 

4.55  e9 

3.95  e9 

-13.2 

Sxy 

Shear  strength  (in-plane) 

7.14  e7 

5.99  e7 

-16.1 

Gxz 

Transverse  shear  stiffness 

8.93  e9 

6.30  e9 

-29.5 

SX7, 

Transverse  shear  strength 

5.76  e8 

4.92  e8 

-14.5 

aThe  mechanical  test  was  performed  by  P.  Moy  of  ARL,  and  a  separate  report  will  be  published. 
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Per  figure  5  and  table  1,  note  that  the  x-direction  and  j-direction  geometries  are  not  completely 
balanced  in  terms  of  aligned  fiber  content.  Thus,  the  properties  are  slightly  different  between 
these  two  directions.  This  is  also  accounted  for  in  the  model.  However,  as  seen  from 
comparison  of  experimental  and  predicted  values,  the  degree  of  variation  is  somewhat 
underpredicted.  In  most  cases,  the  model  and  experimental  values  show  good  agreement.  In¬ 
plane  axial  and  shear  properties  show  an  overall  average  of  90.2%  prediction  accuracy  across 
strength  and  stiffness  predictions. 

Figure  7  shows  the  x-direction  microstress  state  in  the  RVE  under  unit  extension.  As  expected, 
the  tows  aligned  in  the  loading  direction  take  the  bulk  of  the  loading.  The  transverse  stresses  that 
this  load  case  places  on  the  weft  tows  perpendicular  to  the  load  direction  often  lead  to  pull-apart 
or  matrix  microcracking  of  the  tows.  However,  ultimate  failure  is  governed  by  the  strength  of 
the  warp  tows,  which  still  holds  stress  after  the  weft  tows  begin  to  fail. 
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Figure  7.  Contour  plot  illustrating  axial  stress  under  unit  normal  strain  along  the  x-direction 
axis. 


Shear  stresses  in  the  RVE  under  unit  in-plane  shear  strain  are  shown  in  figure  8.  For  this  planar 
loading,  stresses  are  much  more  evenly  distributed  across  the  warp  and  weft  tows.  The  bulk 
matrix  itself  is  also  a  significant  load  carrier,  as  the  inter-tow  shear  stresses  themselves  are  often 
governed  by  the  matrix  properties  (similar  to  the  matrix-driven  shear  properties  of  a 
unidirectional  laminated  composite). 
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Figure  8.  Section-cut  contour  plot  illustrating  internal  in-plane  shear  stress  under  unit  shear  strain.  (Note  that 
two  stacked  RVEs  are  shown.) 


Transverse  shear  stiffness  is  not  well  predicted,  although  the  model  does  predict  transverse  shear 
strength  reasonably  well.  This  is  a  subject  of  current  work.  Even  though  background  derivations 
have  not  been  developed  herein,  it  should  be  noted  that  appropriate  boundary  conditions  to 
represent  transverse  shear  are  not  as  readily  obtainable  as  those  for  other  deformation  states. 
Thus,  the  disagreement  of  model  and  experiment  for  this  property  is  not  unexpected  since  some 
adjustment  to  the  boundary  conditions  corresponding  to  transverse  shear  may  be  needed. 

Nevertheless,  strength  can  still  be  well  predicted.  As  for  a  given  microstress  field,  the  model  can 
adequately  determine  the  maximum  allowables  for  failure  initiation.  In  short,  transverse  shear  is 
not  a  periodic  condition  because  a  state  of  constant  transverse  shear  is  always  accompanied  by  a 
state  of  nonuniform  moment  distribution.  Thus,  if  a  continuum  comprised  of  RVEs  is 
assembled,  each  with  a  uniform  transverse  shear,  the  moment  present  in  each  RVE  across  the 
continuum  will  vary  with  length.  Periodicity  is  disrupted,  and  no  single  displacement  boundary 
condition  can  represent  what  is  seen  in  each  RVE.  Consequently,  the  periodic  boundary 
conditions  representing  the  in-plane  shear  and  axial  strains  (also  applicable  to  flexure  states, 
though  not  treated  in  the  current  work)  may  not  be  similarly  derived  for  the  transverse  shear 
state. 

With  the  ability  to  quickly  evaluate  multiple  stress  states,  the  micromechanics  model  has  been 
used  to  predict  the  in-plane  failure  envelope,  as  shown  in  figure  9.  This  provides  a  useful  “at-a- 
glance”  evaluation  of  the  load  capacity  of  a  given  material  under  a  variety  of  applied  loads.  Also 
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Figure  9.  Plane-stress  failure  envelope  with  comparison  to  common  failure  theories. 


shown  in  figure  9  is  a  comparison  to  the  predictions  of  the  well-known  maximum  stress  and 
Tsai-Wu  failure  theories.  All  failure  theories  are  inherently  formed  to  agree  at  the  uniaxial 
failure  points  (essentially,  strength  properties).  For  this  relatively  simple  loading,  the  general 
agreement  between  the  Tsai-Wu  failure  ellipse  and  the  current  micromechanics  model 
predictions  is  fairly  close,  with  most  disagreement  occurring  under  the  biaxial  loading  of 
quadrants  I  and  III.  (Note  that  Tsai-Wu  or  other  common  failure  theories  do  not  include 
provision  for  applied  transverse  loads  or  bending  moments  in  their  failure  prediction,  whereas 
the  micromechanics  herein  can  accommodate  such  considerations.) 

Considerable  scatter  is  seen  in  the  results  of  the  off-axis  testing — 22%  and  30%  standard 
deviation  for  the  30°  and  45°  specimens,  respectively.  This  is  largely  due  to  the  inherent 
inhomogeneity  of  any  textile  fiber  structure,  especially  with  the  relatively  sparse  through-thickness 
stitch  pattern  of  the  3-D  orthogonal  composite.  Eight  tests  were  conducted  for  each  of  the  30°  and 
45°  specimens  to  help  offset  this  drift,  lend  more  veracity  to  the  average  test  result,  and  maintain 
validity  to  comparisons  with  model  predictions.  For  the  30°  specimens,  the  micromechanics- 
based  model  overpredicted  the  maximum  allowable  peak  load  to  within  6.4%  of  the  average  of 
the  experimental  results.  With  considerably  more  scatter  and  two  suspect  data  points,  the  peak 
load  of  the  45°  off-axis  specimens  was  overpredicted  by  the  model  by  35%.  However,  if  two 
unreasonably  low  suspect  data  points  were  removed,  prediction  of  the  model  agreed  within  16% 
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of  the  experimental  data.  Given  this  testing,  the  micromechanics  model  allowed  effective 
prediction  of  failure  for  this  material  in  simple  and  relatively  more  complex  load  cases. 

3.2  Effects  of  Stitch  Density  on  Stiffness  and  Strength  Properties 

Figures  10-14  illustrate  the  effect  of  stitch  density  on  the  axial  and  transverse  shear  stiffness  and 
strength  of  the  3-D  orthogonal  composite.  Discrete  points  represent  the  model  predictions,  and  a 
solid  line  indicates  a  curve  fit  that  has  been  applied  to  each  plot.  An  R  value  is  shown  on  each 
plot  to  indicate  the  relative  closeness  of  each  curve  fit.  Stitch  density  of  1.6  stitches/cm 
represents  the  initial  RVE,  as  shown  in  figure  1 ,  which  is  then  modified  for  each  new  density  per 
the  specifications  of  table  5. 


The  graphs  in  figures  10  and  1 1  present  the  loss  of  in-plane  axial  stiffness  and  strength  for 
increasing  stitch  density  as  predicted  via  modeling  (axial  properties  in  the  y-direction  are  not 
plotted  because  results  compare  nearly  identically  to  the  x-axis  properties  and  thus  do  not 
provide  new  information).  Due  to  the  increased  presence  of  interruption  of  the  axial  fiber  tows 
which  present  an  in-plane  stress  concentration  and  loss  of  in-plane  fiber  volume  fraction,  it  is 
seen  that  stiffness  and  strength  decrease  linearly  as  stitch  density  is  increased.  Linear  trends  for 
both  property  losses  follow  reasonably  well  with  a  “rule  of  mixtures”  relationship  or  a  volume 
averaging  of  constituent  properties.  The  results  for  zero-stitch  density  have  been  directly 
modeled  within  the  FEM  but  also  compare  closely  to  values  obtained  for  a  cross-ply  laminate 
analysis  in  the  Classical  Laminate  Theory  ( 9 ),  as  the  3-D  weave  essentially  reduces  to  a  cross-ply 
laminate  once  the  z-direction  tows  are  removed. 
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Figure  11.  Model  predicted  effect  of  stitch  density  on  axial  strength. 


stitch  density  (stitch  /  cm) 


Figure  12.  Model  predicted  effect  of  stitch  density  on  through-thickness  (z-axis)  stiffness. 
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Figure  13.  Model  predicted  effect  of  stitch  density  on  transverse  shear  stiffness. 


Figure  14.  Model  predicted  effect  of  stitch  density  on  transverse  shear  strength. 
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Figure  12  plots  the  effect  of  stitch  density  on  through-thickness  stiffness  (z-direction)  of  the  3-D 
weave.  In  typical  applications,  stresses  in  the  thickness  direction  may  not  be  important,  even  for 
thick  composite  parts,  but  this  is  a  property  which  is  directly  affected  by  stitch  density  and  can 
represent  an  important  characteristic  of  structural  capability.  An  increase  in  stitch  density  has  a 
direct  effect  towards  increase  in  through-thickness  stiffness.  As  with  the  axial  properties,  this 
linear  trend  follows  reasonably  well  with  what  would  be  predicted  via  a  rule  of  mixtures  property 
approximation.  A  similar  plot  for  through-thickness  strength  is  not  shown  in  the  current  study 
due  to  a  present  lack  of  experimental  verification  of  properties  and  failure  modes  in  this  loading 
condition.  An  ongoing  side  study  has  shown  that  stitch  pullout  and  tow  peeling  are  the  failure 
modes  in  this  case.  Experimental  methods  for  obtaining  reliable  properties  are  currently  under 
development. 

Figures  13  and  14  show  the  effect  of  stitch  density  on  transverse  shear  properties.  Over  the 
selected  range  of  stitch  densities,  transverse  shear  properties  are  the  most  highly  affected  by 
stitch  density  and  can  offer  90%  or  greater  gains  in  stiffness  and  strength  compared  to  an 
unstitched  sample.  This  is  accompanied  by  a  relatively  lower  magnitude  of  loss  of  in-plane 
properties  on  the  order  of  20%,  clearly  illustrating  the  reasoning  behind  the  microgeometry  of 
the  3-D  orthogonal  weave  and  showing  the  value  of  its  employment  in  impact  loading  or  in  thick 
components. 

Although  it  has  been  shown  that  the  current  micromechanics  model  predicts  transverse  shear 
strength  well  but  does  not  predict  transverse  shear  stiffness  with  complete  accuracy,  the  model 
can  be  calibrated  to  fit  the  experimental  data  and  still  accurately  predict  trends  that  arise  from 
modification  of  the  RVE  microgeometry  for  the  parametric  study.  Also  note  that  if  greater  stitch 
densities  were  investigated  beyond  those  currently  studied,  the  improvement  shown  in  figure  13 
should  approach  a  region  of  diminishing  returns  as  stitch  densities  approach  volume  fractions  on 
the  order  of  the  warp  and  weft  tows. 

The  mechanical  response  under  transverse  shear  loading  is  relatively  more  intricate  than  for  axial 
or  planar  shear,  as  evidenced  by  the  more  complex  deformation  and  stress  profiles  seen  under 
these  conditions.  Consequently,  prediction  of  these  properties  is  not  amenable  to  simple 
approximations  and  must  be  modeled  directly.  Changes  in  transverse  shear  properties  for  altered 
stitch  density  do  not  follow  the  rule  of  mixtures  approximations.  Furthermore,  properties  of  the 
unstitched  specimens  cannot  be  compared  to  predictions  from  the  Classical  Laminate  Theory 
because  as  these  methods  do  not  include  provisions  for  calculation  of  transverse  shear  properties. 
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4.  Conclusions 


In  this  report,  micromechanical  methods  are  used  to  analyze  the  failure  of  a  3-D  orthogonal 
woven  textile  composite.  Photomicrographs  of  the  weave  under  investigation  have  allowed  the 
creation  of  an  accurate  RVE  within  the  FEM.  The  micromechanical  techniques  have  been 
employed  to  determine  the  stiffness  and  strength  of  the  3-D  weave  under  axial,  in-plane  shear, 
and  transverse  shear  loading.  In-plane  axial  and  shear  properties  show  an  overall  average  of 
90.2%  prediction  accuracy  across  strength  and  stiffness  predictions.  Transverse  shear  stiffness  is 
not  as  well  predicted,  although  the  model  does  predict  transverse  shear  strength  reasonably  well, 
with  86%  accuracy  compared  to  other  experiments.  Extending  these  results,  a  parametric  study 
has  been  performed  to  determine  the  effect  of  variation  in  stitch  density  upon  the  consequent 
mechanical  properties.  This  parametric  study,  which  evaluates  the  effects  of  stitch  density  on 
mechanical  performance,  illustrates  the  potential  use  of  the  micromechanics  methods  as  an 
optimization  tool  as  well  as  a  tool  for  material  investigation  and  development.  In-plane  stiffness 
and  strength  are  somewhat  adversely  affected  by  increased  stitch  density,  whereas  through¬ 
thickness  axial  properties  show  improvement.  These  linear  trends  could  be  reasonably  well 
estimated  by  a  rule  of  mixtures  volume-averaging  approximation.  However,  transverse  shear 
properties  exhibit  a  more  complex  mechanical  response,  which  does  not  follow  simple  analytical 
approximation  and  must  be  modeled  directly.  Over  the  selected  range  of  stitch  densities, 
transverse  shear  properties  are  the  most  highly  affected  by  stitch  density  and  can  offer  90%  or 
greater  gains  in  stiffness  and  strength  compared  to  an  unstitched  sample  and  a  20%  relatively 
lower  magnitude  of  loss  of  in-plane  properties. 
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